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PACS 87 . 15 . km - Protein-protein interactions 

PACS 05 . 65 . +b - Self-organized systems 

PACS 87. 16. A — Theory, modeling, and simulations 

Abstract. - We discuss a realistic scenario, accounting for the existence of sub-micrometric 
protein domains in cell membranes. At the biological level, such membrane domains have been 
shown to be specialized, in order to perform a determined biological task, in the sense that they 
gather one or a few protein species out of the hundreds of different ones that a cell membrane 
may contain. By analyzing the balance between mixing entropy and protein affinities, we propose 
that such protein sorting in distinct domains can be explained without appealing to pre-existing 
lipidic micro-phase separations, as in the lipid raft scenario. We show that the proposed scenario 
is compatible with known physical interactions between membrane proteins, even if thousands of 
different species coexist. 
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Introduction. — Membrane functional organization 
js a ubiquitous issue in cell biophysics It has be- 

come consensual that membrane constituents, predomi- 
nantly lipids and proteins, adopt a non-random, hetero- 
geneous organization |3j. Lateral segregation is now ac- 
cepted as a fundamental requirement for membrane bio- 
logical functions [4jj7| . The immense variety of lipids and 
proteins in a single biomembrane (several hundred differ- 
ent species) leads to a large variety of interactions be- 
tween them. These interactions have been demonstrated 
to favor the formation of membrane domains, whose size 
ranges from few nanometers to microns. Domains can be 
induced by lipid-lipid interactions (from which the concept 
of "raft" emerged [2|4]|8]), lipid-protein interactions [9 -12 
or protein-protein interactions [13f|17| . Understanding the 
role of these domains and their physico-chemical origin re- 
mains a key problem in cell biology. 

The physical scenario investigated in this work belongs 
to the protein-protein category above (even though lipids 
do play a role because they indirectly participate in ef- 
fective inter-protein interactions 10 - 12] 18 19]). The pro- 
posal that protein-protein interactions can drive the for- 
mation of domains independently of a lipidic micro-phase 
separation has recently been advanced by several research 
Statistical mechanics arguments have 



while short-range attraction favors condensation of mem- 
brane proteins in a dense phase, some weaker repulsion at 
longer-range, the origin of which is still debated, prevents 
a complete phase separation. The resulting phase at equi- 
librium is called a "cluster phase" 15 , by analogy with 
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similar phases in soft condensed matter 

However, these studies consider the clusterization of a 
single species of proteins, without taking into account the 
hundreds of different species in the cell membrane with 
which it coexists. In particular, they do not explain why 
these other species are excluded from the clusters under 
consideration. A priori, mixing entropy would stabilize 
multiple-species clusters (Fig. [T]). Understanding the spe- 
cialization of clusters remains a matter in debate of cell 
biology. Experiments demonstrate the existence of finite- 
size protein domains gathering a single or few species. 
They are believed to have a precise biological function, 
e.g. (high-fidelity) signaling [5 13 21 - 24 , cell adhesion 



and motility |25| , immune response 
exocytosis and membrane fusion 
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26 , endocytosis 27 
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groups 

been proposed, relying on the same global mechanism: 



The present work proposes to solve this issue by cou- 
pling the previous cluster-phase model 15 16 - account- 
ing for a finite domain size - and a Flory-Huggins the- 
ory 29 - taking into account the great variety of mem- 
brane proteins. Proteins are classified into families so that 
the same-family proteins, which are not necessarily iden- 
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Fig. 1: Two schematized views of membrane cluster phases, 
with q — 4 different protein families (different colors). Top: 
if proteins of different families have a sufficiently high affinity, 
they mix together in clusters. Bottom: in the opposite case, 
clusters are mono-colored, i.e. they are highly concentrated in 
proteins of one family because the ensuing energy gain is larger 
than the related entropic cost. This phase-separation leads to 
cluster specialization because proteins that are destined to co- 
localize in same clusters indeed do so, spontaneously. 



tical, are those that tend to co-localize in same domains 
to form functional platforms, e.g. a G-protein-coupled re- 
ceptor and associated G proteins and effectors [5||6 22 
The references 



13 14 17 21 23-28 provide other exam- 



ples. The model assumes that the same-family proteins 
have a higher "affinity" at contact than proteins in dis- 
tinct families, i.e. contact between proteins in the same 
family is energetically favored. This difference of affinity 
is measured by a parameter x > 0, reminiscent of a Flory 
interaction parameter. 

We prove that if x is smaller than a critical value x c ; 
clusters mix protein families because entropy dominates. 
If X > X C i clusters demix and are essentially composed of 
a single family, the different families being segregated be- 
tween different specialized clusters (Fig. [I]) . Furthermore, 



X c is small even if the number q of different families is 
very large (up to 10 4 ) because x c ^ ksTlnq, where fc^T 
is the thermal energy. The main result of the present work 
is that a contact energy difference as small as 2 to 3 kgT 
suffices to favor a demixed cluster phase where clusters are 
specialized because they gather specific proteins. We show 
that such contact-energy differences are indeed compatible 
with known interactions between membrane proteins. 

Description of the model. We consider N 
isotropic (or weakly anisotropic) proteins interacting in 
a 2D medium that accounts for the lipid "sea" in which 
proteins reside. Even though they are not explicitly taken 
into account, lipids are responsible for effective interac- 
tions between proteins, of elastic or entropic nature (see 
below). We assume that the total interaction free energy 
of a protein configuration is the sum of pairwise potentials, 



L/jv(ri,...,rjv) 



l<i<j<N 



(1) 



where is the position of protein i. Un does not depend 
on the protein orientations because we assume that they 
are isotropic or that anisotropy is sufficiently weak (see 
Refs. [30j|3l] for the quantification of anisotropy effects) ; 
Uij depends on i and j, because proteins can be of different 
nature^] Assuming a pairwise interaction is certainly an 
approximation since many-body effects are known to exist 
in this context 30 , 32 . 33 , 36 . Thus the present calculation 
is only a first step towards the full solution. We shall 
return to this point in the Discussion section. 

Following previous works (see 14 15 17 
therein), we assume that the potentials 
tive at short range (R < 1 run 



and references 
u are attrac- 
i.e. roughly speaking 
"at contact") and weakly repulsive at intermediate range 
(R > 10 nm). The repulsion has been given several expla- 
nations Mj 15j 17 34 which we do not intend to discuss 
here (see Discussion). Indeed, in the "liquid-droplet" for- 
malism adopted below, the precise shape of the potential 
does not qualitatively affect the results, in agreement with 
numerical work |15| . Note that even a weak (a fraction of 
fcfiT) intermediate-range protein-protein repulsion suffices 
to promote clusters [i~5] . 

As compared to these previous studies, we introduce 
a new ingredient in our model: the dependence of the 
short-range part of Uij on protein families i and j. 
Origins of the energy modulation at contact are man- 
ifold: direct electrostatic and polar interactions, even 
though screened in physiological conditions, play a role at 
short range; hydrophobic-mismatch interactions depend 
on transmembrane-protein hydrophobic thicknesses 35- 
|37| ; in membrane with different lipid species, proteins 
recruit lipids in their immediate vicinity for which they 
have a higher affinity, because of electrostatic and hydro- 
gen-bond interactions between lipids and polar or charged 
amino acids of proteins. These lipid annuli, which "wet" 



Note that in this paper, all energies and thermodynamic poten- 
tials are implicitly expressed in units of kgT. 
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the proteins, are also responsible for effective attractive 



forces depending on the protein nature 10 12 



Results. — We first focus on the case where Uij is in- 
dependent of i and j (identical particles). We recall briefly 
the formalism and notations introduced in Ref. [16] (see 
also [38]). The canonical partition function for identical 
indistinguishable particles is 

A-iiW f 

Z i N ) = ^rT dr 1 ...dr Jv e- t ^ ri >- r ^, (2) 

iV! JyN 

where the length A making Z dimensionless is set to be 
the particle diameter UM. V is the d-dimensional vol- 
ume where particles evolve. Speaking of membrane nano- 
domains assumes that one is able to define regions V k 
of V that partition the N particles into disjoint clusters: 
there are N\ monomers, N 2 dimers, and so forth, so that 
N = kNk ■ Each fc-particle cluster (or fc-mer) dwells in 
a distinct region V k . Then neglecting interactions between 
the different clusters 34 38 , the integral in Eq. ([2| can be 
written as a product of integrals over the V k . After simple 



algebraic manipulations 16 , 



{N k } k 



-F(k) 



(3) 



where we have introduced the free energy of a fc-cluster, 
F(k): F(i) = for monomers and for fc > 1, 



F(k) 



In 



A <i(i-fc) 
fc! 



dri . . .dr fe _i e 



-u k 



-fk + -fk 3 + a k a . 



(4) 
(5) 



The last approximation, valid when k 3> 1, is at the core 
of the analytical calculations in Ref. [16] and will also be 
adopted in the present work [^] Equation ^ is a gener- 
alization of the liquid-droplet model, to which a repulsion 
term is added (see Ref. [16] for further details) : (i) the first 
term accounts for the short-range attraction between pro- 
teins and — / < is a bulk free energy per particle; (ii) the 
second one represents the free-energy cost of the free in- 
terface between the cluster and the surrounding fluid and 
7 > is a line tension. In dimension d — 2, this term scales 
as Vk, which means that protein domains are assumed to 
have a disk shape [15], Other shapes, such as stripes, can 
exist at very high protein concentrations that are not con- 
sidered in the present work [39]; (iii) the last term, where 
a > 1, takes into account the weak longer-range repul- 
sion, the strength of which is measured by a > 0. Since it 
dominates at large fc, it renders too large clusters unsta- 
ble and is responsible for their finite size at equilibrium. 



2 Note that as compared to this expression, k was replaced by 
— 1 in Ref. [16] for technical reasons, which has no consequences 
in the large k limit of our interest, because it amounts to corrections 
of the order of 1/k. The present formulation will be more easily 
generalized below. 



Note that the exact form of F for small multimers is irrel- 
evant because they are virtually nonexistent in the regime 



of parameters at hand 16 

In practice, we are interested in cases where cluster 
sizes range from a few dozen to a few thousand parti- 
cles. This implies that a < 1 (weak intermediate-range 
repulsion) and /, 7 3> 1 (moderate short-range attrac- 
tion) [16] . For simplicity, we also set a = 3/2, and 
d = 2 because a membrane is two-dimensional. How- 
ever the present derivations could be extended to any 
d and a. The mean volume fraction of fc-clusters, c&, 
derives from Z(N) after introducing the chemical po- 
tential /i and switching to the grand-canonical ensem- 
ble: c k = (N k )A d /V = exp[^fc - F(k)] = exp[-G(fc; /x)], 
where G(k; fx) = F(k)-fj,k = -(/ + /z)fc + 7fc 1/2 + crfc 3 / 2 is 
the grand potential of a fc-cluster [16] . The expected value 
of the protein volume fraction, 4> = (N)A d /V, sets the 
chemical potential /i through the condition <p = £\ =1 kc k - 
At low volume fraction (ft, fi = In c\ < In <j> takes large neg- 
ative values and G(k) increases monotonously with fc: c k 
is maximal at fc = 1 and decreases exponentially with fc, 
with typical width . The system is essentially com- 
posed of monomers. Above a critical concentration <f> c , the 
chemical potential \i > [f is such that G(k) has two local 
minima: one at fc = 1 and one at fc* 3> 1. The proteins 
are partitioned between a gas of monomers and clusters 
of typical aggregation number fc*. At the critical poten- 
tial /i c , G(fc;/i c ) has an inflexion point at fc = fc* = fc c . 
Thus //= and fc c satisfy d k G(k c ; /i c ) = d 2 k G(k c - n c ) = 0. 
Consequently [f — — f + y/3aj and k c — 1 + 7/ {3a) ^> 1. 

In the general case, where the system contains q families 
of proteins, one has to introduce new definitions. Mx is 
the number of proteins of family K and N = ^ Mk ■ 
If two particles of families K and K' are adjacent in the 
same cluster, their binding energy at contact is denoted by 
ck,k' ■ Due to thermal agitation, an entropic contribution 
must be added to this energy, leading to an average free 
energy per bond lpk,k> ~ £k,k> 16 . This free energy can 
be embodied in the Flory interaction parameter Xk,K' = 
v<Pk,k' < [29], where v is the average number of particle 
neighbors in a cluster bulk (v = 6 in the present case of 
2D dense clusters). If xk = M K /N, the free energy per 



particle of a homogeneous mixture of q families reads 29 



F 



K ~ Yl X K^X K 
K=l 



E Xk,K'X K x k > 

K,K' = 1 



(6) 



Before tackling the general case with q families of pro- 
teins, we first focus on q — 2, where A- and _B-familics co- 
exist. To simplify the discussion, we assume that xaa = 
Xbb and define x = Xab — Xaa > 0. The canonical par- 
tition function becomes Z(Ma,Mb) and Eq. ^ must be 
adapted to this new situation: if N kA kB is the number of 
clusters containing kj± A-proteins and fcg B-proteins, then 

Z= E II {J ^-r r 1 • (7) 



,} kA,kE 



k A ,k t 
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Here 0) — F(0, 1) = for monomers and 



F(k Al k B ) = - In 



A d(i-fe) 
M fa! 



fe-i 
k i=i 



i e 



(8) 



if fc = ItA+kB 3> 1- Simplifying this expression to a liquid- 
droplet-like form [Eq. ([5])] requires to discuss further the 
internal organization of clusters. Two main scenarios must 
be envisaged [40] : either A- and -B-proteins are demixed in 
a same /c-cluster, with an interface between both phases, 
or they are homogeneously mixed. In contrast to previ- 
ous studies where clusters are tackled in a canonical for- 
malism because they cannot exchange material with their 
environment 40 , the demixed case is always unfavorable 
in the present case for the following reason: it implies an 
energetic cost proportional to y/k per cluster due to the 
interface between A- and B-rich phases. Thus replacing 
an assembly of demixed clusters by an assembly of A- and 
B-rich ones, of mixed type, is always favorable, because 
the ensuing entropic cost is of order 1 per cluster, thus 
<C Vk . Therefore we only consider mixed clusters be- 
low. If one defines the ^4-protcin fraction of a /c-cluster, 
x = kA/k, the Stirling formula leads to 

F(k A ,k B ) = k[xlnx + (l-x)\n(l-x)} (9) 
■ A d(i-k) 



In 



fc! 



dri . . .dr fc _i e" 



for k 3> 1. The first term is the mixing entropy. The 
second one can now be written in a liquid-droplet-like form 
and Eq. Q becomes 

F(k A , k B ) = -fk + -fk 1 ' 2 + afc 3 / 2 + kh 2 (x) (10) 

where — / = xaa/2 and 

h 2 (x) = x^(l — x) + xlnx + (1 — x) ln(l — x). (11) 

This expression of h 2 is only a mean-field approximation 
but it will become exact in the large q limit of interest 
below 41 . Furthermore, we have chosen a line tension 7 
independent of x and this is also an approximation. We 
shall return to this point below. We now introduce two 
chemical potentials ha and hb and G(kA,kB\ Ha, Hb) = 
F(kA, ks) — Ha^a — Hb^b- Writing fc^ = xk and ks = 
(1 — x)k, minimizing G with respect to kA and ks amounts 
to minimizing it with respect to k and x. Beginning with 
x, we remark that h 2 (x)— ha x ~ ^b{^~ %) is the mean-field 
free energy of the Ising model in a magnetic field (hb — 
Ha)/^- Here we are interested in the situation where no 
protein family dominates, i.e. where Ma ~ Mb or ha — 
Hb- Thus we focus on the symmetric case ha = Hb = H 
(or Ma = Mb)- Then the mean-field Ising model in a 
vanishing magnetic field presents a second-order transition 
at x° = 2- For a weak difference in protein affinities (x < 
X c ), x x — 1/2 is the most probable value and clusters 
typically contain an equal number of A- and £?-proteins. 
By contrast, if x > X°i h 2 {x) has two non-trivial minima 



at x x and 1 — x x . When x grows beyond Xc clusters 
are essentially mono-colored. Once x x is determined, G 
becomes a function of k: 

G(k; h) = -[/ + M - h 2 (x x )]k + 7 fcV2 + a j,3/2 (12) 

Therefore the problem of minimization of G with respect 
to k now amounts to a one-family problem as considered 
previously, with a renormalized /: f T = f — h 2 (x x ) and 
thus a renormalized /i c . Note that if x = 0, then x x = 1/2 
and / r = / + ln 2. Indeed in the two-color case, the canoni- 
cal partition function Z is normalized by Ma ! Mb ! instead 
of N\ because same-color particles are indistinguishable. 
Thus the entropy is increased by In 2 per particle as com- 
pared to the one-color case (when Ma = Mb)- 

The generalization to any q is straightforward: F and 
G become functions of k and of the number fractions 
xk = kn/k (K — l,...,q). We focus again on the 
most symmetric case: xk,k — Xo is independent of K, 
Xk,k> — Xm if K 7^ K' and M K = N/q. We define 
X = Xm - Xo > and 



h q (x u x 2l . . . ,x q ) = x ^2 X K%K' 

K<K' 



^ x K \nx K , (13) 



K=l 



the mean-field free energy of the g-state Potts model [29] , 



exact in the large q limit 41 
when q > 3. On gets 



The transition is first-order 



X c = 2^ ln(g-l)~21n 9 . 
q- 2 



(14) 



At x c ; the majority number fraction is x c — 1 — 1/q, 
other colors have identical fractions x' — l/[q(q — 1)], and 
h c q ~ — 1/q. Again, clusters are essentially mono-colored 
at the transition and beyond, while colors are mixed if 
X < x c - Figure [2] shows the typical phase diagram derived 
from Eq. ( |14[ ), and from the fact that is renormalized: 
Hi- = fi +~3aj, with / r = f — h q . Note that if x = 0, 
/ r = / + lnq for the same reason as above. It must be 
emphasized that in this phase diagram the "transitions" as 
determined above are not true thermodynamic transitions 
because the finiteness of clusters smoothes transitions 142] . 
More technically, the system grand potential J (still in 
units of k^T) can be calculated: 



J 
V 



X! C ki,-,k K 
k 1 ,...,k K 



(15) 



where c kl ,...,k K = exp[-G(fc x , . . . , k q ; Hi, ■ ■ ■ ■ ,U q )] is the 
mean volume fraction of (fci, . . . , fc g )-clusters |^] The fact 
that G is dominated at large k by ak a with a > en- 
sures the uniform convergence of the series and of all its 
derivatives, and thus prevents any singularity. 



3 Eq. ([l5j| implies that PV = -J = N cl k B T, where P is the d- 
dimensional pressure and N c \ = V^2 c k lt ...,k K is the total number 
of clusters. Thus the system behaves like an ideal gaz of clusters. 
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Fig. 2: Phase diagram in the (fi, x) plane. It is calculated here 
in the case q — 3 but its qualitative aspect is valid for any q, 
with Xc — 21ng [Eq. (14 1]. The lines represent crossovers and 



not true thermodynamic transitions, as discussed in the text. 

Note however that the real critical value differs from 
above when q > k, because the xk do not take contin- 
uous but discrete values, xk — kn/k in a fc-cluster (kx 
an integer). Thus no more than k variables xk can be 
simultaneously non-vanishing and h q is to be minimized 
on this subset of R d , which amounts to a minimization 
of hk- Thus \ c — 21nfc and x c ~ 1 — 1/k. Generally 
speaking, if k* still denotes the typical cluster size, then 
X c — 21n[min(fc*, q)]. In a similar way, in the phase dia- 
gram of Fig. [2j q should be replaced by min(k*,q). 

Discussion. — By appealing only to generic two-body 
protein-protein interactions, we have shown that protein 
nano-domains may spontaneously form in biomembranes 
at equilibrium. One of the goals of the present approach is 
to account for the finite size of membrane nano-domains. 
Out-of-equilibrium arguments have also been developed to 
explain this finite size 43 44 . In this respect, our goal here 



has not been to discuss the relative merits of equilibrium 
or out-of-equilibrium approaches, which are certainly com- 
plementary, but to show that the equilibrium approach 
alone can yield finite-size nano-domains. Exploring the 
additional effects of out-of-equilibrium phenomena in the 
cluster phase scenario could clarify their role. In partic- 
ular, taking into account membrane recycling 43 44 or 



proteins interacting differently with their partners when 
in an active state 



formalism. 

Furthermore, we have shown that, in spite of the associ- 
ated entropic cost, it becomes favorable to demix different 
protein families in distinct domains when x > X c with 
Xc — 21n[min(fc*,q)]. Here k* is the typical cluster size. 
Even if min(fc*,g) ~ 10 3 to 10 4 , a realistic upper-bound 
for a real plasma membrane, x° ~ 14 to 18, which means 



that \e K ,K - e K ,K'\ - \<PK,K ~ ¥K,K<\ = X°/ v ~ 2 to 
3ksT. Is it realistic to expect such differences of con- 
tact energy between different protein families? The pre- 
cise calculation of interaction energies for proteins embed- 
ded in a lipidic membrane is a tedious task that can be 
hardly tackled analytically. By contrast, several numer- 
ical studies have estimated the variations of the energy 
at contact for transmembrane proteins with variable hy- 
drophobic thicknesses 18 19 36 37 . All these studies are 



consistent with contact energy differences larger than 3 
ksT. As far as lipid-mediated interactions are concerned, 
we are not aware of any quantitative results in the liter- 
ature. However, given the order of magnitude of this in- 
teraction |11||12| , its amplitude is likely to be modified by 
one ksT or more due to interplay between lipid and pro- 
tein species. Thus typical modulations of \tx,K — e K,K'\ 
can easily exceed 2 or 3 fc^T. At the biological level, 
this result shows that an important concentration of pro- 
teins in a biological membrane leads to their gathering and 
sorting in nano-domains. By contrast, below the critical 
concentration, they would be distributed randomly on the 
membrane. To this respect, increasing the concentration 
strongly facilitates the encounter of proteins having a con- 
tact energy slightly lower than the average one, and thus 
favors biological functions. 

The line tension 7 measures the free energy cost for a 
protein to be at the cluster boundary as compared to the 
bulk. Indeed, a bulk particle has about twice as many 
neighbors as one at the boundary. Thus one would expect 
7 to depend on x x and thus on x- But clusters are essen- 
tially mono-colored at the transition and beyond, and one 
expects 7 to depend only weakly on x m this region of the 
phase diagram. 

To derive the form of the cluster free-energy F(kA, ks) 
in Eq. (10 1, or its g-color counterpart [see Eq. (fl3|], we 
have assumed pairwise interaction potentials in Eq. ([I]), 
because inferring F^kAjks) in this particular case is eas- 
Real potentials contain many-body contributions 



ler 



1G 



30 33 36 



that have for example been explored in Refs 
In another work, long-range repulsion has been proposed 
to arise from steric interactions, of many-body origin by 
nature 14 . Even though it is a complex issue, fully taking 



45 should be tractable in the present trix 33 



into account many-body interactions will be necessary to 
derive properly the effective parameters /, 7, a or x from 
first principles. Numerical simulations, that are out of 
the scope of the present work, will certainly be required, 
because, e.g., the full calculation of long-range iV-body 
potentials mediated by the elastic membrane requires to 
compute the inverse and the determinant of a N by N ma- 
The estimation of many-body effects associated 



with hydrophobic mismatch also requires intensive simual- 
Note however that our whole argument relies 



36 



tions 

on the effective free energy F(k 1 , . . . , k q ) and not on the 
exact shape of the potential C/jy from which F(k\ 1 . . . , k q ) 
ensues. 

It would be interesting to investigate the role of strongly 
anisotropic inclusions in the future, because in this case, 
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the long-range interaction can become attractive 31 . We 



anticipate that the cluster-phase scenario remains valid if 
only a minority of proteins are strongly anisotropic, that 
are homogeneously mixed in clusters (we recall that pro- 
teins of a same family are not necessarily identical but are 
classified according to their short-range affinities ; thus 
they are miscible in each cluster). Indeed, in this case, 
the repulsive energy still grows faster than the cluster ag- 
gregation number k (i.e. a > 1 in Eq. ([5])), because the 
majority of long-range interactions in the cluster remain 
repulsive. The main ingredient for the existence of stable 



clusters is preserved 15 



We have also assumed that \K,K' can only take two 
possible values, \o if K = K' °r Xm if K ^ K'. The real- 
ity is more complex because of the great variety of mem- 
brane proteins. Future studies would also have to con- 
sider more realistic distributions of ^-parameters around 
these typical values. The asymmetric case where Mk (or 
Hk) depends on K should also be explored. The function 
G{{kn }; {fiK }) will then have to be minimized numeri- 
cally, providing distributions of the xk and cluster sizes k. 
However, we anticipate that our conclusions will be qual- 
itatively unchanged since the present approach captures 
the essential physical ingredient, namely the competition 
between energy and configurational entropy. 

Finally, our predictions ought to be tested experimen- 
tally. Fdrster Resonance Energy Transfer (FRET) is an 
appropriate tool to quantify the distance between tagged 
bio-molecules 



4G 



After incorporating different pro- 
teins species in giant vesicles, tagged with different flu- 
orophores [22], it would be possible to test qualitatively 
the predictions of the present work by playing on physical 
parameters such as the asymmetry of the proteins or the 
thickness of the lipidic membrane. 



I am indebted to Manoel Manghi and Laurence Salome 
for helpful discussions and comments, and to Revaz Ra- 
mazashvili for his kind reading of the manuscript. 

REFERENCES 

[1] Lenne P.-F. and Nicolas A., Soft Matter, 5 (2009) 2841. 
[2] Jacobson K., Mouritsen O.G. and Anderson R.G., 

Nat. Cell Biol, 9 (2007) 7. 
[3] Engelman D.M., Nature, 438 (2005) 578. 
[4] Simons K. and Ikonen E., Nature, 387 (1997) 569. 
[5] Davare M.A., et al., Science, 293 (2001) 98. 
[6] Laporte S.A., Oakley R.H. and Caron M.G., Science, 

293 (2001) 62. 

[7] Rozenfeld R. and Devi L.A., Trends Pharmacol. Sri., 

31 (2010) 124. 
[8] Pike L.J., J. Lipid Res., 47 (2006) 1597. 
[9] Sprong H., van der Sluijs P. and van Meer G., Nat. 

Rev. Mol. Cell Biol, 2 (2001) 504. 
[10] Poveda J. A., et AL., Biochim. Biophys. Acta, 1778 

(2008) 1583. 



[11] Gil T., Sabra M.C., Ipsen J.H. and Mouritsen O.G., 

Biophys. J., 73 (1997) 1728. 
[12] Gil T., et al., Biochim. Biophys. Acta, 1376 (1998) 245. 
[13] Daumas F., et al., Biophys. J., 84 (2003) 356. 
[14] Sieber J.J., et al., Science, 317 (2007) 1072. 
[15] Destainville N., Phys. Rev. E, 77 (2008) 011905. 
[16] Destainville N. and Foret L., Phys. Rev. E, 77 (2008) 

051403. 

[17] Gurry T., Kahramanogullari O. and Endres R.G., 

PLoS One, 4 (2009) e6148. 
[18] de Meyer F.J.M., Venturoli M. and Smit B, Biophys. 

J., 95 (2008) 1851. 
[19] West B., Brown F.L.H. and Schmid F., Biophys. J., 

96 (2009) 101. 
[20] Stradner A., et al., Nature, 432 (2004) 492. 
[21] Park P.S. and Palczewski K., Nat. Chem. Biol, 1 

(2005) 184. 

[22] Abankwa D. and Vogel H., J. Cell Sri., 120 (2007) 
2953. 

[23] Prior LA., et al., Nat. Cell Biol, 3 (2001) 368. 
[24] Tian T., et al., Nat. Cell Biol, 9 (2007) 905. 
[25] Espenel C, et al., J. Cell Biol, 182 (2008) 765. 
[26] de Bakker B.I., et al., ChemPhysChem, 8 (2007) 1473. 
[27] Frick M., et al., Curr. Biol, 17 (2007) 1151. 
[28] Hess S.T., et al., Proc. Natl. Acad. Sci. USA, 104 
(2008) 17370. 

[29] de Gennes P.-G., Scaling concepts in polymer physics 

(Cornell University Press, Ithaca, 1979). 
[30] Dommersnes P.G. and Fournier J.-B., Eur. Phys. J. 

B, 12 (1999) 9. 

[31] Chou, T., Kim K.S., and Oster G., Biophys. J., 80 
(2001) 1075. 

[32] Kim K.S., Neu J.C., and Oster G.F., Europhys. Lett, 
48 (1999) 99. 

[33] Fournier J.-B., Dommersnes P.G. and Galatola P., 

C. R. Biologies, 326 (2003) 467. 

[34] McConnell H.M. and Moy V.T., J. Phys. Chem., 92 
(1988) 4520. 

[35] Dan N., Pincus P. and Safran S.A., Langmuir, 9 
(1993) 2768. 

[36] Brannigan G. and Brown F.L.H., Biophys. J., 92 
(2007) 864. 

[37] Schmidt U., Guigas G., and Weiss M., Phys. Rev. 
Lett, 101 (2008) 128104; Phys. Rev. Lett, 102 (2009) 
219802. 

[38] Mitchell D.J. and Ninham B.W., J. Chem. Soc. Trans. 

2, 77 (1981) 601. 
[39] Sear R.P., et al., Phys. Rev. E, 59 (1999) R6255. 
[40] Christensen A., Stoltze P. and Norskov J.K., J. 

Phys.: Condens. Matter, 7 (1995) 1047. 
[41] Mittag L. and Stephen M.J., J. Phys. A: Math. Nucl. 

Gen., 7 (1974) L109. 
[42] Ruckenstein E. and Nagarajan R., J. Phys. Chem., 

79 (1975) 2622. 
[43] Foret L., Europhys. Lett., 71 (2005) 508. 
[44] Turner M.S., Sens P. and Socci N.D., Phys. Rev. 

Lett, 95 (2005) 168301. 
[45] Phillips R., Ursell T., Wiggins P. and Sens P., Na- 
ture, 459 (2009) 379. 
[46] LAKOWICZ JR., Principles of Fluorescence Spectroscopy, 

3rd edition (Springer- Verlag, New York, 2006). 



p-6 



